Long-Range Magnetic Interactions in the Multiferroic Antiferromagnet M11WO4 
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The spin- wave excitations of the multiferroic MnWCU have been measured in its low-temperature 
collinear commensurate phase using high-resolution inelastic neutron scattering. These excitations 
can be well described by a Heisenberg model with competing long-range exchange interactions and a 
single-ion anisotropy term. The magnetic interactions are strongly frustrated within the zigzag spin 
chain along c axis and between chains along the a axis, while the coupling between spin along the b 
axis is much weaker. The balance of these interactions results in the noncollinear incommensurate 
spin structure associated with the magnetoelectric effect, and the perturbation of the magnetic 
interactions leads to the observed rich phase diagrams of the chemically-doped materials. This 
delicate balance can also be tuned by the application of external electric or magnetic fields to 
achieve magnetoelectric control of this type of materials. 

PACS numbers: 78.70.Nx, 75.30.Ds, 75.85.+t 



Magnetoelectric multiferroic materials, which exhibit 
the coexistence of ferroelectric (FE) and magnetic orders, 
have attracted great attention in recent years i"— Sev- 
eral classes of multiferroics among transition metal oxides 
have been discovered including geometrically-frustrated 
CuM0 2 (M is Fe, Cr)££ RbFe(Mo0 4 ) 2 7 , Ni 3 V 2 OA 
or rare-earth (R) manganites i?Mn03 and i?Mn 2 05*2r— 
The ability to simultaneously control the electric (E) 
and magnetic (M) properties makes these multiferroics 
promising candidates for technological applications^^ 
A characteristic feature in those magnetically induced 
multiferroics is the presence of long-range magnetic struc- 
tures with noncollinear spiral spin configurations. Such 
magnetic order is a consequence of magnetic frustration 
either due to geometric constraints or competing ex- 
change interactions resulting in a close competition of 
different magnetic structures that are nearly degenerate 
in energy. 

The mineral Hiibnerite MnWCU appears to be an 
unique material that not only exhibits intriguing mul- 
tiferroic phenomena but also shows rich magnetic phases 
via chemical substitutions J^~— It has been considered 
one of the prototypical multiferroics capable of magne- 
toelectric (ME) controlJ^ Unlike i?Mn0 3 where the 
spiral magnetic structure often involves ordering of the 
rare-earth moments, MnWCU is a frustrated antiferro- 
magnet (AF) with only one type of magnetic ion. The 
Mn 2+ spins (5=5/2) undergo successive transitions in 
zero field* 2 ^ The low-temperature (T) magnetic structure 
has a collinear spin configuration [Fig. 1(a)]. For T be- 
tween 7.8 K (T/vi) and 12 K (T/y 2 ), the magnetic struc- 
ture evolves into an incommensurate (ICM) elliptical spi- 
ral configuration accompanied by a spontaneous electric 
polarization P along the crystalline b axis. When T is 
further raised between T N2 and T N3 (^ 13.5 K), MnW0 4 
becomes a collinear ICM and paraelectric. 



In MnW04, the electric polarization that is correlated 
with the spiral magnetic structure can be well under- 
stood by the microscopic picture regarded as inverse 
Dzyaloshinskii-Moriya interaction^ - — On the other 
hand, characterizing the magnetic interactions that cause 
the formation of the complex spin structures and un- 
derstanding how the modification of exchange couplings 
affects the evolution between different phases remains 
an unresolved issue despite intensive experimental and 
theoretical studies j^ 7 - - — Early inelastic neutron scatter- 
ing (INS) work suggested that the stabilization of the 
collinear configuration requires higher order magnetic 
interactions. 27 Although the magnetic exchange param- 
eters obtained by Ehrenberg et al fit the experimental 
data, the longest bond distance was associated with the 
strongest coupling constant. Later, density functional 
calculation and classical spin analysis were performed to 
investigate the magnetic structure and FE polarization 
in MnW04^ Tian et al. concluded that the spin ex- 
change interactions are frustrated along both the a and 
c axes. However, a quantitative experimental character- 
ization of the magnetic interactions is still lacking. Here 
we report high-resolution INS measurements that show 
that the low-T magnetic ground state of pure MnWC^ 
indeed results from the competition of long-range inter- 
actions that are highly frustrated and sensitive to small 
perturbations. The comprehensive mapping of the of the 
magnetic excitations along several symmetry directions 
allowed an unambiguous determination of the dispersion 
relations and the exchange interactions. Most impor- 
tantly, such microscopic characterization of the spin cou- 
pling constants provides a fundamental step toward the 
construction of the ground state Hamiltonian from which 
the FE phase can be derived. 

A 5-g single crystal of MnWCU was grown by the 
floating-zone technique. Neutron diffraction was per- 
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FIG. 1: (Color online) (a) The magnetic structure of MnWCU 
in the collinear, commensurate phase at low temperature. 
The magnetic spins lie in the ac plane with the moment 
canted to the a axis about 35°.— The magnetic spins form 
zigzag fN4- chains along the c axis, and are coupled antiferro- 
magnetically along the b axis, (b) The magnetic interactions 
along and between spin chains in the be plane, (c) Higher 
order magnetic interactions along the a axis direction. The 
magnetic couplings are labeled with increasing bonding dis- 
tance. Note the monoclinic crystal structure (/3 = 91.14°) 
makes J6/J7, J10/J11 pairs different. 



formed on a small piece of this crystal (0.2 g) to ver- 
ify the spin structure using the four-circle single crystal 
diffractometer HB3A at the High Flux Isotope Reactor at 
the Oak Ridge National Laboratory (ORNL). The inelas- 
tic neutron scattering measurements were performed us- 
ing the Cold Neutron Chopper Spectrometer (CNCS) at 
the Spallation Neutron Source at ORNL. The momentum 
transfer wavevectors q = (q x ,Q y ,Qz) are in units of A -1 
at positions (H,K,L) = (q x a / 2it , q y b / 2it , q z c / 2it) in re- 
ciprocal lattice units (rlu), where a — 4.83 A, b = 5.75 A, 
c = 4.99 A. We aligned the crystal in several scat- 
tering planes such that the spin- wave (SW) dispersion 
along the [1,0,-2], [1,0,2] and [0,1,0] symmetric directions 
that pass across the magnetic Bragg peaks can be read- 
ily measured. The incident neutrons with wavelength of 
A = 4.4 A were chosen to ensure the needed energy reso- 
lution to separate various magnetic excitation branches. 

MnW04 orders in the collinear phase with two inequiv- 
alent commensurate wavevectors qM = (1/4, 1/2, ±1/2). 
Figures 2(a) and 2(b) present the spin excitation spec- 
tra along the [1,0,2] direction with K = 0.5 and 1.0. 
This scanning direction goes through the magnetic Bragg 
peak (1/4, 1/2, 1/2). The data clearly show four distinct 
branches that disperse out from the magnetic zone cen- 
ter (ZC) to the zone boundary. The spectra reveal a spin 
gap of 0.5 meV and boundary energy around 2.2 meV. 
The excitation bandwidth is consistent with the energy 





FIG. 2: (Color online) (a) SW dispersion spectra along the 
[1,0,2] direction through the magnetic peak (1/4,1/2,1/2). (b) 
Magnetic excitation spectra along the same direction with 
K — 1.0. (c and d) Calculated spectra along [1,0,2] direction 
with K — 0.5 and K — 1.0 using magnetic exchange interac- 
tions up to Jn with instrument resolution convoluted. The 
solid lines overlapped with experimental data in (a) and (b) 
are the predicted dispersion curves. 



scale of the ordering temperature of 13.5 K. The inten- 
sity of the excitation spectra is highly asymmetric with 
respect to the magnetic Bragg point. For example, the 
spectral weight of lowest branch in Fig. 2(a) is completely 
missing as H approaches zero, while it shows the highest 
intensity as it moves toward H = 0.5. This highlights 
the importance of a complete survey in reciprocal space 
to fully map out the magnetic dynamics. 

Figures 3(a) and 3(b) display the magnetic scattering 
spectra along the [1,0,-2] direction that crosses the other 
magnetic wavevector (1/4,1/2,-1/2), while Figs. 4(a) 
and 4(b) illustrate the spectra along the [0,1,0] direc- 
tion. The icattering intensity map shows similar asym- 
metric feature on both sides of the magnetic ZC. Investi- 
gations along those two directions have been previously 
reported^ but the lack of spectral weight for certain 
branches makes a correct description of the SW disper- 
sion difficult. 

The SW dispersion curves can be modeled by a general 
effective Heisenberg Hamiltonian: 
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where . indicates summation over pairs of spins f 
D is the single-ion anisotropy, and S{ Z denote the spin 
components along the easy axis. To calculate the corre- 
sponding spectral weight, the low-T spin structure has 
been verified by collecting a complete set of magnetic 
reflections that covers the entire reciprocal space. The 
spin configuration obtained by Rietveld refinement with 
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TABLE I: Magnetic exchange coupling parameters derived from spin wave model calculation according to the Eq. 1. The Mn 2+ 
ions located at position (1/2, 0.685, 1/4) interact with neighboring spins through one or two oxygens. The bonding distance (in 
units of A) between Mn. . .Mn are also listed. The magnetic interaction constants from previous experimental and theoretical 
studies have been normalized in units of meV for comparison. 
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FIG. 3: (Color online) (a) SW dispersion spectra along the 
[1,0,-2] direction through the magnetic peak (1/4,1/2,-1/2). 
(b) Magnetic excitation spectra along the same direction with 
K — 1.0. (c)-(d) Calculated spectra along the same symmet- 
ric directions with K — 0.5 and K — 1.0. 



FIG. 4: (Color online) (a) SW dispersion spectra along the 
[0,1,0] direction through the magnetic peak (1/4,1/2,-1/2). 
(b) Magnetic excitation spectra along the same direction 
through the nuclear peak (0,1,0). (c)-(d) Calculated spec- 
tra along the same symmetric directions with H — 0.25 and 
H = 0. 



group theory analysis is in good agreement with a previ- 
ous report. 23 This configuration is then used to evaluate 
the magnetic scattering cross-section: 

^ cx f(Q)e- 2W ~ Q a Qp)S ap (Q,co) (2) 

where f 2 (Q) is the magnetic form factor for the Mn 2+ 
ion, e~ 2W is the Debye- Waller factor, Q a is the a compo- 
nent of a unit vector in the direction of Q, and S a P(Q, uo) 
is the response function that describes the aj3 spin-spin 
correlations. 

Combining data along all symmetry directions, the dis- 
persion relations can be simultaneously modeled using 
Eq. (1). We started with nine exchange parameters and 
single-ion anisotropy D, as assumed in Ref. 27. Although 
the fitting parameters provide a fair description of the 
data along the [1,0,-2] and [0,1,0] direction, they fail to 
capture the dispersion relations along the [1,0,2] direc- 



tion. Instead, an extra pair of magnetic interactions Jio 
and Jn have to be included to achieve satisfactory agree- 
ment for the data along all scan directions. As listed in 
Table I, the magnetic exchange constants generally de- 
crease in amplitude as the bonding distances increase, 
except the weaker J2 and J5 (both close to zero) along 
the b axis. Mn 2+ has an electronic configuration of 3d 5 
(orbital singlet) that is not expected to have any mag- 
netic anisotropy, and the nonvanishing D = 0.09 meV in- 
dicates a possible spin-orbit coupling that causes the pin- 
ning of the magnetic moments in the ac planed To test 
whether those parameters are consistent with the actual 
spin structure, the magnetic energy for all possible spin 
configurations (2 8 = 256 with eight spins in one magnetic 
unit cell) are calculated. We verified that only the spin 
structure depicted in Fig. 1(a) gives the lowest energy. 
In addition, we employed the spin-rotation technique to 
calculate the expected intensity for the SW modes 31 and 
to compare with the observed wavevector dependence of 
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the spectra weight. The right panels of Figs. 2-4 show the 
excitation spectra map using Eq. (2) with fitted parame- 
ters. The excellent agreement between the calculated in- 
tensities and the experimental data provides convincing 
evidence that the formation of the collinear spin struc- 
ture indeed requires longer range magnetic interactions. 

A complex spin configuration would form in the con- 
ventional ID frustrated spin chain if the next-nearest- 
neighbor antiferromagnetic (AFM) interaction J' be- 
comes substantially stronger compared to the nearest- 
neighbor (NN) ferromagnetic (FM) exchange coupling 
J. A spiral phase might appear for J' / \J\ > 1/4 and a 
collinear ttil structure will emerge when J f / \J\ > 1/2^ 
The zigzag E-type phase observed in i?Mn03 is a clas- 
sic example of competing short-range FM-AFM inter- 
actions caused by lattice distortions^^ In the case 
the MnW04, the exchange coupling are predominantly 
AFM and three-dimensional. Nevertheless, strong long- 
range magnetic interactions comparable to NN interac- 
tion (Ji) along the c axis (\Ja/ J\\ > 1/2) and a axis 
(l^sMI > 1/2, \Je/Ji\ ~ 1 and \J 7 /Ji\ > 1/4) are ob- 
served. This reflects the magnetic frustration in those 
two directions and is consistent with the ICM compo- 
nents present in the a and c directions but not in the 
b direction, when the system enters the spiral phase. 
The exchange coupling remains sizable (Jio) even at a 
rather long distance. Such an unusual extended inter- 
action could also be viewed as a much-reduced NN ex- 
change coupling because of a nearly 90° Mn-O-Mn bond- 
ing ang le^ Thus, the unique crystalline structure makes 
MnW04 a promising material to achieve novel physical 



properties when the magnetic interactions are modified. 
It was reported that doping a few percent magnetic or 
nonmagnetic impurities can drastically affect the spin 
order iZr— For instance, while replacing Mn 2+ with Fe 2+ 
ions that have a larger local magnetic anisotropy seems 
to enhance the collinear structure^ the introduction of 
nonmagnetic Zn 2+ ions that weaken the overall magnetic 
interactions switches the ground state from collinear to 
spiral order j 18 i 2Q Those results demonstrate that chemi- 
cal substitutions are a viable tool to tune the multiferroic 
properties in the extremely sensitive MnWC^. 

In summary, high-resolution INS is used to study the 
SWs in the collinear phase of MnWC^. The collinear 
spin order is stabilized by the delicate balance of com- 
peting long-range magnetic interactions. We provide an 
effective Hamiltonian to describe the highly frustrated 
magnetic order within the zigzag spin chains along the c 
axis and between spin chains along a axis. Rich and com- 
plex magnetic phases are expected in chemically doped 
MnW04 due to the fine-tuning of the magnetic inter- 
actions. The delicate balance of exchange interactions 
in MnW04 can also be used to achieve magnetoelectric 
control using external electric or magnetic fields. 
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